Interferometric distillation and determination of unknown two-qubit entanglement 
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We propose a scheme for both distilhng and quantifying entanglement, applicable to individual 
copies of an arbitrary unknown two-qubit state. It is realized in a usual two-qubit interferometry 
with local filtering. Proper filtering operation for the maximal distillation of the state is achieved, 
by erasing single-qubit interference, and then the concurrence of the state is determined directly 
from the visibilities of two-qubit interference. We compare the scheme with full state tomography. 
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Introduction. — Multiparticle interference is a strik- 
ing phenomenon connecting with quantum entanglement. 
For pure states, the connection is rather straightforward. 
In a two-particle interferometry_[l|, 0| , the interference 
visibility gives the concurrence [3, |j| j 8' widely used en- 
tanglement measure, of the two-particle pure state [5j . In 
a multiparticle Aharonov-Bohm interferometry [yj, the 
visibility can be used to prove the quantum nonlocal- 
ity of Greenberger-Horne-Zeilinger entanglement For 
mixed states, however, multiparticle interference comes 
from a mixture of entanglement and classical correlation, 
and it is hard to distinguish the two different correlations. 
It is interesting to find a way to extract entanglement 
from the interference, which is the aim of this work. 

In quantum information research, there are strong de- 
mands of distilling and quantifying entanglement Q. 
Currently available schemes are of two types, one us- 
ing multiple copies of a target state and the other using 
individual copies. Since the multiple copies are harder 
to prepare in laboratory in general, it may be necessary 
to explore further the latter type. The distillation of 
the latter type has been done using local filtering, for 
a known two-qubit state [l^l or after full state tomog- 
raphy And no scheme of the latter type has been 
proposed for directly quantifying an entanglement mea- 
sure, such as concurrence, of an arbitrary mixed state 
in experiments; note that the existing schemes of the for- 
mer type for determining concurrence have not been real- 
ized or provide a lower bound of concurrence for 
mixed states, while concurrence was recently determined 
in experiments by using two copies of a pure state [T3 |. 
Therefore, it is valuable to find a scheme of the latter 
type for distilling and directly determining entanglement 
of an unknown state (without full state reconstruction). 

In this work, we propose an interferometric scheme for 
both distilling and determining entanglement, applicable 
to individual copies of an arbitrary unknown two-qubit 
state. It can be realized in a two-qubit interferometry 
with local filtering [10, iH. The maximal distillation (the 
normal form ilSj) of the state is first achieved, by itera- 
tively erasing single-qubit interference, and then the con- 
currences of both the initial and the distilled states are 
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FIG. 1. Two-qubit interferometry with local filtering. It con- 
sists of a source, distillation parts, and detection parts. The 
source generates individual copies of an arbitrary unknown 
state of two qubits A and B, each represented by pseudospins 
tj and ij, J = A, B. Qubit A files to the detectors (Da) of 
Alice, passing through three beam splitters (BS), while B to 
Bob (Db); see solid arrows. The state is transformed into its 
maximally distilled state in the distillation parts, and then 
its concurrence is determined by measuring the visibilities of 
two-qubit interference in the detection parts. 



directly determined from the visibilities of two-qubit in- 
terference. This quantification is based on our important 
findings that the two-qubit interference shows three dif- 
ferent "local" extrema (visibilities) in general and that 
when the single-qubit interference is fully erased, the 
three extrema give the Lorentz singular values [15|, a 
linear combination of which gives the concurrences. Our 
scheme is conceptually different from full state tomogra- 
phy and practically useful. 

Two-qubit interferometry with local filtering. — We in- 
troduce the interferometry (Fig. [1]). The source gener- 
ates individual copies of a state of two qubits A and B, 
which fly to Alice and Bob, respectively. The qubit de- 
gree of freedom, the pseudospin, can be photon polariza- 
tion, particle path, particle spin, etc. For illustration, we 
choose particle path as the pseudospin, by considering 
two particles A and B, each injected to either the upper 
(pseudospin up) or the lower path (down) of its side. 

The 4x4 density matrix p of the initial two-qubit 
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state is represented by using pseudospin basis states, and 
written by using a real matrix R as |16j] 

/5=^ ^ RwCFi®av, Rw =TT{pai®av), (1) 

where Tr(- • ■ ) means the trace of matrix, (Tq is the 2x2 
identity matrix, and (T/'s {I = 1,2,3) are the Pauli ma- 
trices. Then the single-qubit states {pj = Tijp) of Alice 
and Bob are represented by 4- vectors {Rqo, Rio,R20, R30) 
and (i?oO) ^01) ^02, ^03), respectively; Trj(---) means 
the trace over the degrees of freedom of qubit j j). 

In the distillation parts, which are absent in usual 
interferometries [ll, 0, Alice (Bob) has local operation 

Dj{fj,0dis.j,4'dis,]) = F{fj)U{0dis,j,4>dis,j) on qubit j = 
A (B). It transforms the initial state p into its normal 
form [ll,[i3| Pdis, its maximally distilled state p^ . 

Pdis = (£'a®£'b)/5(^a®£'b)^ (2) 

The local filtering F{fj) and the rotation tJ{0dis,j,4'dis.,j) 
of qubit j are supported by two beam splitters, BS2J and 
BSij, respectively, and represented as 

^(/) ^ (J ?) , ^(M, ^ -if ) ^ 

(3) 

Here, < /j < 1 is the filtering parameter controlled by 
the reflection amplitude of BS2J, ^disj G [Oi""] parame- 
terizes the transmission at BSij, and ^disj £ [0j27r] is 
the phase shift. In the filtering operation F{fj), parti- 
cle j is abandoned with probability 1 — /J when it flies 
along the lower path after scattering by BSij. Whether 
qubit i is filtered off or not is certified at detector D3J. 
The two beam splitters constitute the minimal setup for 
the distillation. This is understood from the fact [l^ 
that each local operation on qubit j corresponds to a 
Lorentz transformation of the 4- vector of qubit j . F and 
U correspond to the Lorentz boost and the spatial rota- 
tion, respectively. We emphasize that the Lorentz boost 
mathematically introduced in Ref. [15] is physically re- 
alized here by F, using beam splitter BS2j. We will see 
later how Dj is efficiently found for an unknown state p. 

In the detection parts, Alice and Bob count the num- 
ber Tiij of the particles j arriving at detector Dij, i ~ 1,2, 
during measurement time, long enough to get the statis- 
tical average of coincidence correlation {nip^nii-o) for a 
given setting of all the beam splitters. They tune BS3A 
and BS3B to see single- and two-qubit interferences in 
(nij(nij-t-rt2j)) and (5niA(5niB) = (niA"iB)-(niA(«iB+ 
"■2b))('T'1b(?^ia + '^2a)), respectively; the other correla- 
tions involving D2j- contain the same information. Here, 
the number riij is normalized by the total number N of 
states ending in neither D3A nor D3B, and ('T-ij"(f^ij+'T-2j)) 
means the average number of qubit j in the events 
where the other qubit j (7^ j) is not filtered (does not 
end in D3J). The qubit rotation at BS3J is represented by 

U{9dct,j-,4>dct,]), where 6'dctj G [0,7r] and (j)dct,j e [0, 27r]. 



Note that the phase accumulation of particle j along its 
path is absorbed in the rotation angles 4>dis,j and 4>dct,j- 
The visibilities of single- and two-qubit interferences 
are defined, respectively, as [T9| 

Vj=A,B = W[{nxj{n^j + n2j))], 
Vab = VK[(,5niA<5niB) + l/4], 

where W[x\ = (max[a;] — min[a;])/(max[x] -I- min[a:]) 
and max[x] (min[x]) means the local maxima (min- 
ima) of X over the parameter space of BS3A and BS3B, 
{6'dct,A,'/>dot,A,6'dct,B,</>dot,B}- Since the mean value of 
((5niA'5niB) over the space is zero, the ad hoc factor 1/4 
is added in Eq. so that the values of Vab are equal to 
the local maxima of {SnixSniB) ■ Note that one needs to 
tune the transmission probability cos^ -^j^ of BS^w, to 
obtain the information of the diagonal parts of p [19| . 

Entanglement distillation. — We first explain how to 
transform p into pdis- It is based on the facts [l3| that 
local density matrices pdis,j = Trjpdis of the normal form 
are proportional to the identity matrix and that any ini- 
tial state p can be transformed iteratively to its normal 
form by filtering operations. In our interferometry, we 
find that these facts are realized as follows. When pdis is 
achieved in the distillation parts, the single-qubit inter- 
ference visibilities Va and Vb vanish, since the identity 
generates no interference. Thus, one can obtain Dj in an 
iterative way such that Alice and Bob alternately tune 
her/his beam splitters of the distillation parts to make 
the visibility Vj of her/his single-qubit interference van- 
ish, until Va and Vb both vanish simultaneously. 

We describe each iteration step. In the (2fc— l)-th step, 
k — 1, 2, • • • , Alice observes the single-qubit interference 
signal in (niA(«iB + "^2b)) by tuning BS3A, with setting 
her parameters as (/A,^^dis,A) — (I7O) but without tun- 
ing (/B^'°"^^^d?s'',B^^'^^Sf',B^^) fixed by Bob; in the first 
step, Bob starts with (/B°\^d'isB) — (I'O). By compar- 
ing the signal with its general form, (niA('T'iB + ?12b)) = 
(1 -I- 7A • wa)/2, Alice determines 7a = (7ai,7A2,7A3), 
which in fact represents the spatial part of the 4-vector 
of qubit A; the general form has only a pair of ex- 
trema ±|7a|, i-e., Va (= |7a|) is single-valued. Here 
Vj = (sin 6'dctj- cos (t)dct,j , sin 6'dctj sin 0dct j , cos 6'dctj ) is 
the rotation vector of BS3j. Then, by setting f'^'' = 
\/ {1 ~ |7a|)/(1 + |7a|) and the rotation vector of BSia in 
the (2A; — l)-th step as w^^^^"^^ = — 7a/|7a|, Alice achieves 

the situation that V^^*" vanishes; after the setting, 7a 
rotates to be parallel to z= (0, 0, —1) at BSia, and then 
vanishes by the filtering at BS2A- Next, Bob performs his 
(2fc)-th step in the same way as Alice's {2k — l)-th step, 
except for the exchange A -f-^ B. After such iterations, 
the distillation parameters converge to (/j, ^disj, 0disj), 
at which Va = Vb = and pdis is obtained. 

Entanglement determination. — Before discussing en- 
tanglement quantification, we first show an interesting 
feature of Vab- For a state p', transformed from p by 
an arbitrary setting of the distillation parts, we derive 
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a compact form of the cross-correlation, (SniASniB) = 
jV^Q'v^, where Vj is the rotation vector of BSsj, the 
colmnn vector t/g is the transpose of vb, Q' is the 3x3 
matrix defined by Q[, = i?;,,/Ko ^ Vrn'/W^ 
1,1' — 1,2,3, and R is the real parametrization of p' 
in Eq. ([T]); the number normalization by N gives the 
factors l/i?oo and l/(i?Qo)^ in QJ^,. From this compact 
form and the fact that Vj^A.B spans over the surface of 
unit sphere, it is easy to see that ((^tiia^'^ib) has three 
pairs of "local" extrema ztAj's (Ai > A2 > A3 > 0), i.e., 
Vab has the three values A;'s, and that Aj's are identical 
to the singular values of Q' up to sign factor. Here, Ai 
is the global maximum of ((5niA<5niB), while A2 (A3) is 
the maximum over the space of va and vb orthogonal to 
VA,i and vb.i {va,i, vb,i, va,2, and '!Tb,2), where Vj^i is the 
rotation vector of BS3J at which (SniASniB) shows A;. 

The above finding becomes very useful when pdis is 
achieved in the distillation parts (p' — pdis and R — 
Rdis). In this case, the visibilities A/'s give the Lorentz 
singular values of the initial state as 

N 

■So = , r J J ' Sl=SoAl, S2 = S0A2, S3 = gsoA3, (5) 
Ai Ja/b 

since they are equal to the singular values of Rdis / (/a /b ) , 
Rdis,io = -Rdis.oi = 0, and -Rdis, 00 = N/M. Here, M is 
the total number of injection of p from the source, q = 
Det(i/A,i,WA,2,WA,3)Det(wB,i,WB,2, WB,3) is the sign factor 
guaranteeing the correct singular value decomposition, 
Det(- • • ) means matrix determinant, and {vj^i,Vj^2,'Vj^3) 
is the matrix whose columns are Vj^s. Using the relation 
[Til l between concurrence and Lorentz singular values, we 
find an important result that the concurrences C of p and 
Pdis are directly obtained from Vab , 

C{p) = soC(pdis), 

1 (6) 
C(/5dis) = max[0, -(-1 + Ai + A2 - qh)]- 

Examples. — In Fig. [2l the concurrence is determined 
at each A:-th iteration step, for two examples of p, using 
Vab and Eq. For typical cases of p (non-asymptotic 
case) [Fig. UKa)], vj^^^ g vanishes rapidly within a few 
steps, and the determined value C'^^^p) approaches to the 
exact value C{p) more rapidly. In this case, the deviation 
ofC^'^H/S) fromC(p) is estimated 20] as \C{p) ~C'^''\p)\ oc 
(vj*^-*)^ for small vj*^-* (< 0.1). Thus, one can determine 
a precise value of C even before the complete distillation. 
The properties of particular types of p are given below. 

(i) When p is pure, only one distillation step is neces- 
sary, since V^'^' = v'^'' for all k due to the complementar- 
ity [5^. Note that {vfj^l^f+C^ip) = 1 for /a = /b = 1- 

(ii) When p is separable and uncorrelated, p = pa®Pb, 
the local properties of A and B are independent. There- 
fore, only two steps are required, Vab = 0, and C{p) = 0. 
Particularly, when either pA or pe is pure, its single- 
qubit visibility is one at fc = 0, and p cannot be distilled 
as Pdis vanishes. When p is separable but has classical 
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FIG. 2. (Color online) The concurrence C*'^'(/3) (dot-dashed 
lines) , determined from the two-qubit visibilities Vab , and the 
single-qubit visibilities V^'°' and V^*"' in the fc-th distillation 
step are shown for (a) /3,,a = \\(i>e){(j)e\ + ^(\ ti)(ti I + 
1 itXit I) with (e,A) = (0.5,0.8) and (b) p = \\ tt)(n 
+ + where = (e| ft) + I ii 

))/^/m2 and |*±) = ^(| it) ± | ti))- The case (a) is 
typical, showing rapid convergence to p^ia, while (b) is an 
asymptotic case with slow convergence. For comparison, the 
exact value i4] of C(/5) is given (dotted lines). 



correlations, on the other hand, more than two steps are 
neccessary, Vab 7^ 0, and C(p) =0. 

iii) When p is a Werner state or a Bell-diagonal state 
|2l|, the distillation is not necessary, 
(iv) There is the so-called asymptotic case [l^, where 
large number of steps are necessary [Fig. EJb)] and most 
states are abandoned by the filtering _s. g). 

Below, we propose an optimal way of determining con- 
currence. After the distillation, one first determines all 
the rotation vectors Vj^i, by observing (^niA^J^iB) in nine 
difi^erent settings of BS3A and BS3B and comparing the 
results with the compact form of {SriiASriiB) (derived be- 
fore). Then, one measures ((5niA<5niB) (thus A;) at the 
determined setting of BS3j=A(B) at va(b),i- We empha- 
size that a crude determination of Vjj is enough for a 
precise detection of A; and C. It is because A/ is a local 
maximum, around which a small error (^ 5) in the direc- 
tion Vj^i for A/ causes only a much smaller error 5"^) in 
the value of A;. This makes our scheme efficient. Table I 
shows that for states not much filtered, our scheme is as 
efficient as the tomography for the quantification of 
the initial state. For the distillation and the quantifica- 
tion together, it can be more efficient than previous tomo- 
graphic schemes, e.g., in Ref. [lj[; the previous schemes 
require roughly 2-4 times larger number of state copies 
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State Distillation Quantification Total Tomography 

Bell 2400 (0) 9 x 100 + 3 x 200 3900 360 

Werner 2400 (0) 9 x 500 + 3 x 4000 18900 38700 

1/4 2400 (0) 9 X 100 + 3 X 7800 26700 24300 

Po.e 2400 (0) 9 X 300 + 3 X 3800 16500 30600 

po.9,0.6 7200 (1) 9 X 500 + 3 X 6800 32100 54900 

po.5,0.8 40800 (5) 9 X 600 + 3 X 9600 75000 38700 

TABLE I. Monte Carlo simulation of the minimum number 
of individual copies of a given state p, that need to be used 
to determine its concurrence [or so(— 1 + Ai + A2 — qX3)/2 in 
Eq. (O] within ±0.01 statistical error in our scheme (fourth 
column) and by full state tomography [S^l (fifth). In our 
scheme, it is the sum of the number of necessary copies for 
the distillation (second column) with fcdis iterative steps, and 
that for the quantification (third) consisting of nine measure- 
ment settings for the determination of vj^i and three set- 
tings for the three maxima A;'s; the parentheses show fcdis- 
In the distillation, the copies are used to achieve and to 
check vj'''^'°' < 0.1; for the tested states, this condition of 
yj'=dis) jg enough to obtain C within the ±0.01 error. On the 
other hand, among available tomography schemes, we con- 
sider here the most efhcient one with nine measurement set- 
tings and four detectors. We test six representative states 
usually tested in entanglement detection [a, [13]) ^ Bell state 
1*0) = ^(1 n)±l U)), a Werner state = ||*o){*o|±||, 

1/4, p; = |*o)(*o| ±2^(1 nXU I + I U}{n \), and (in- 
troduced in Fig. [2]) with two different sets of (e, A). Note that 
the efficiency of our scheme strongly depends on p (as the to- 
mography) and becomes worse for states more filtered (those 
with smaller N/M); N/M = 1 (no distiUation; fcdis = 0) for 
the first four states, 0.95 for po.g.o.e, and 0.42 for po. 5,0.8- 



than our scheme, as they require the tomography twice 
(once before and once after the distillation). Moreover, 
our scheme improves previous distillations [13, [HI , as it 
is applicable to unknown states. Therefore, our scheme 
is practically useful, in the situation ^ i23(| that for un- 
known states, the existing schemes are less efficient than 
the tomography and virtually require it. 

Conclusion. — We have proposed a "quantum entan- 
glement concentrator" , in which the entanglement of an 
arbitrary unknown two-qubit state is distilled and deter- 
mined. We remark the following meaningful features. 

First, our scheme is within experimental reach and ap- 
plicable to generic types of qubits, as it has only local 
operations using a tunable beam splitter, currently avail- 
able 10]. Second, we show that even for mixed states, 
concurrence and Lorentz singular values are directly and 
experimentally accessible, interestingly from the extrema 
of two-qubit interference; concurrence has been deter- 
mined experimentally only for a pure state [l^ . This 
motivates to study the features of the singular values 
[2^ 1. Third, entanglement quan tification can be closely 
related with distillation [Tol. l25l [26l|. In our scheme, the 
former can be done after the latter. Finally, our scheme 
may be practically useful (e.g., for teleportation jl^l), as 
it achieves the distillation and the quantification within 
one framework. It would be valuable to generalize our 
scheme to larger systems of multiple qubits, where to- 
mography error estimation becomes less feasible. 
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